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NASA’s Marshall Space Flight Center (MSFC), in conjunction with the United States 
Geological Survey (USGS) and aided by personnel from the Asiromaterials Research and 
Exploration Science group at Johnson Space Center (ARES-JSC), is implementing a new 
data acquisition strategy to support the development and evaluation of lunar regolith 
simulants. The first analyses of lunar regolith samples by the simulant group were carried 
out in early 2008 on samples from Apollo 16 core 64001/64002. The results of these analyses 
are combined with data compiled from the literature to generate a reference composition 
and particle size distribution (PSD) for lunar highlands regolith. In this paper we present the 
specifics of particle type composition and PSD for this reference composition. Furthermore, 
we use Figure of Merit (FoM) routines to measure the characteristics of a number of lunar 
regolith simulants against this reference composition. The lunar highlands regolith reference 
composition and the FoM results are presented to guide simulant producers and simulant 
users in their research and development processes. 


I. Introduction 

T HE current lunar architecture calls for the establishment of a permanently manned lunar outpost by 2020. To 
support this effort, the lunar engineering community is developing and testing technologies for In Situ Resource 
Utilization ( ISRU). excavation and drilling, and mitigation of hazards to machinery and human health. The scarcity 
of Apollo samples for testing purposes necessitates the use of Umar regolith simulants as proxy materials. The 
development and evaluation of simulants requires a highly detailed understanding of regolith particle characteristics. 
The Marshall Space Flight Center (MSFC)-led group has developed Figure of Merit (FoM ) algorithms and compiled 
a theoretical reference highlands regolith material against which to quantitatively compare simulants. 
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If. Figures of Merit 

The purpose of the Figure of Merit (FoM) is to quantitatively compare, for a defined parameter or set of 
parameters, a regoiith simulant to a reference material, where the reference may be a regolith sample, a hypothetical 
regolith sample, or other batches or types of simulants 1 '. Though FoM algorithms exist for particle type (modal) 
composition, particle size distribution (PSD), particle shape distribution, and density, we here compare materials 
only using the particle type composition and PSD algorithms. 

The FoM is defined as a value between 0 and 1. As the two materials become more similar the FoM approaches 
!. Mathematically, the composition FoM is the normalized difference of the composition vectors of the two 
materials (reference and simulant) subtracted from unity. Normalization forces the difference of two composition 
vectors to lie between 0 and 1, and subtraction from unity results in a Figure-of-Merit of 1 for a perfect match to 0 
for no match at all. The FoM for composition may be interpreted as the fraction of material that is the same in both 
materials 1 . 

Mathematically, the particle size FoM is the difference between the size relative frequency distributions, which 
are the derivatives of the size cumulative frequency distributions of the two materials, subtracted from unity. Again, 
a Figure* of- Merit of 1 is a perfect match and 0 is no match at all 1 * . 


in. Apollo Regolith Reference Material 

NASA’s current lunar architecture calls for initial outposts to be placed in the lunar polar regions, which are, to 
the best of our knowledge, mantled by highlands-type lunar regolith. Of the Apollo missions, only Apollo 16 visited 
a site located entirely in the lunar highlands, and so we use Apollo 16 regolith samples as the reference materials to 
guide simulant development and evaluation. Cores are our preferred sample types as they provide samples across a 
depth profile, and many proposed operations on the moon will involve excavation of lunar regolith to depths of tens 
of centimeters. As the lunar architecture expands or evolves, we can incorporate reference materials from other lunar 
locales like maria and the KREEP-enriched Procellarium terrane. 

The reference material used here for Figure of Merit calculations is the integrated suite of subsamples from 
Apollo core 64001/64002, which is a sample of lunar highlands regolith from Apollo 16 Station 4 •••■ the 
southeastern-most site on the flanks of Stone Mountain. We chose the 64001/64002 core because it is a complete 
and intact core, it is deemed representative of Apollo 16 site regolith’’ 4 . and it has been well-studied’' 1 ’. 

A. Reference Particle Type Composition 

All particle type composition data can be found in Table 1 - it includes the reference composition derived from 
literature’’ 4 and analyses as described below as well as the simulants described in the next section. 

1. Literature 

We averaged modal particle type data from two studies of core 64001/64002’’ 4 for the basis of the FoM lunar 
reference material: for sample 6400 T’, the lower -30 cm of the core and from 64002 4 , the top -30 cm of the core. 
Each study examined six size fractions front 20 to 500 pm of six suhsantpies of the core at -5 cm interval. The 
authors use a consistent particle classification system' and calculate a weighted average, by weight % of the size 
fractions, of the compositions of each subsample. We combined these subsample averages to a single mean particle 
type composition of the 20-500 pm portion of the 64001/64002 core. 

This particle type classification 7 has primarily been used with data generated by optical microscopy of very fine 
particles and thus some minerals are not classified to the level of specificity we desire. For instance, pyroxenes are 
not differentiated to c!im> or orthopyroxene, and all spinel minerals (chromite, spinel, and ulvospinel), ilmenite, and 
sulfides are undifferentiated as “opaques”. 

More than 90% of the particles by weight of most lunar regoiith samples are less than 500 pm in diameter 8 . An 
average of -20 wt.% of most regolith falls below 20 pnT, but modal data for this fraction are scarce. Therefore, we 
consider this 20-500 pm dataset to be the best available in the literature for our purposes. 

2. Scanning electron microscopy (SEM) and energy dispersive Xray spectroscopy (EDS) 

We generated modal data from QEMSCAN 4 SEM-EDS microbeam analysis of Apollo 16 samples from drive 
core 64001/64002^*“. The analyzed lunar samples were thin sections 64002,6019 (5. 0-8.0 cm depth) and 
64001,6031 (50.0-53. 1 cm depth.) and sieved grain mounts 64002,262 and 64001,374 from depths corresponding to 
the thin sections, respectively. We analyzed four size fractions from each grain mount sample: 500-250 urn, 1 50-90 
pm, 75-45 pm, and <20 pm fractions. These data are not particle type modal data but rather total area modai% by 
phase, either as type of mineral or glass.*' The area modal % is derived from focused beam analyses at stepped 
intervals. 
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For the lunar reference composition, we use the ratios of certain mineral classes from these SEM/EDS data to 
augment our particle type modal data from the literature. For instance, when the literature data 3 '" report only 
“pyroxene”, we subdivide these into clinopyroxene and orthopyroxene based on the Q E M S C A N*-ge» e rated ratio 0 . 
Furthermore, we divide their “opaques” into ilmenite, Fe-sulfide, and spinels (not further differentiated). 

3. Plagiodases composition 

Plagioclase feldspar is the only mineral for which we currently evaluate chemical compositional variability in the 
FoM algorithm. We use the generally accepted composition of An95 i3 for lunar highland regolith plagioclase. This 
means the plagioclase is 95 molar% of the CaAl 2 Si 2 0 8 (anorthite) end-member and only 5 molar% of the NaAlSi 3 0 8 
(albite) end-member. We include plagioclase composition in the FoM because: 

• it is the most abundant mineral in the highlands regolith'; 

• it is the only mineral for which we have reasonable compositional data in both the regolith and the 
simulants; and 

• lunar highlands plagioclase is more calcic than almost any terrestrial plagioclase, and we view closeness to 
lunar plagioclase composition as a significant marker of simulant fidelity. 

B. Reference Particle Size Distribution 

The particle size distribution data (PSD) for 64001/64002 is taken from the literature 8 , it is an average of 12 
subsamples by weight% of each size fraction. Figure 1 shows the average Apollo 3.6 64001/64002 PSD 8 . All data is 
from dry sieving. 



Sieve Size (microns) 


Figure 1. Particle size distribution data for Apollo 16 core 64001/64002. Data from Graf, 1993 s . 


IV. Simulant Compositions and Particle Size Distributions 

As part of the ongoing effort to evaluate and guide lunar regolith simulant development, we have gathered data 
on existing and previously used simulants so as to compare them to our Apollo Reference material. We have 
analyzed simulants of varying types and sources. These include: JSC-1, JSC-1A, and JSC-iAF (fines) are mare 
simulants approximating Apollo 14 regolith 9 - 12 ; Japan’s FJS-I analyzed here is a mare simulant also approximating 
Apollo 14 regolith, but other variations of FJS-I included materials added to simulate Apollo II and Apollo 16 
regolith'"; MLS- 1 (Minnesota Lunar Simulant) as tested here (processed in a plasma stream to generate a glass 
component) is a high-Ti mare simulant that approximates Apollo II htgh-Ti mare regolith; Canada’s OB-1 is a 
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highland simulant intended for mechanical/physical simulation of regolith; and the NASA/USGS medium- and dust- 
sized NU-LHT-1M, NU-LHT-2M, and NU-LHT- ID are highland simulants intended for general purpose. 


A. Simulant Composition 

/, SEM/EDS analyses 

All of the simulants are considerably less texturaliy complex than the analyzed lunar regolith and so we were 
able to obtain consistent particle type analyses by QEMSCAN* SEM/EDS analysis. The software used for textural 
analysis and particle identification is the iDiscover 4.2 package developed by Intellection, Ltd. and incorporated into 
QEMSCAN te technology. The user is able to differentiate and classify basalt fragments in some simulants (grouped 
as “Isthic fragments” in our FoM. analysis), and the pseudo-agglutinate fragments (grouped as “agglutinates” in our 
FoM analysis) in the NU-LHT series simulants, 

2. Piagioclase composition 

Piagioclase compositions in the simulants, either from published or presented analyses 12 or from best estimates 
based on feedstocks' 2 , are included in the composition FoM analysis. 


Table 1 . Particle type modal data, and piagioclase molar% Anorthite, for the lunar reference material 
and regolith simu lants. See text for data so u rces. . 



64001/ 

64002 

NU- 

urr- 

1M 

NU- 

LHT 

-2M 

OB-1 

JSC-1 

JSC- 

IA 

JSC- 
1 A.F 

FJS-I 

MLS-1 

Lithic Fragments 

31.11 

0.00 

0.00 

0.00 

90.92 

90.92 

91.93 

80.18 

52.28 

Glass 

8 . 88 

22.37 

7.17 

52.63 

0.00 

0.00 

0.00 

0.53 

36.57 

Agglutinates 

32.51 

29.02 

23.49 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Piagioclase 

23.32 

38.78 

54.89 

43.95 

1.54 

1.54 

3.39 

14.11 

2.60 

(Plug. An%) 

95 

80 

80 

75 

68 

70 

70 

50? 

47 

Olivine 

0.00 

2.88 

9.5 1 

0.04 

5,63 

5.63 

4.13 

1.13 

0.0 1 

Clinopyroxene 

0.64 

2.04 

3.98 

0.07 

1 .33 

1 .33 

0.42 

1.20 

2.21 

Orthopyroxene 

3.24 

4.37 

0.20 

0.00 

0.0 i 

0.01 

0.0 f 

0.04 

0.03 

Spinel minerals 

0,03 

0.05 

0.01 

0.19 

0.00 

0.04 

0.02 

0.05 

0.03 

Fe-sui fide 

0.0! 

0.00 

0.04 

0.00 


0.00 

0.00 

0.00 

0.00 

Ca-phospates 

0. ! 2 

0.00 

0.43 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Ilmen ile 

0.13 

0.33 

0.19 

0.00 

0.00 

0.08 

0.00 

0.15 

1 07 

Native Iron 

0.01 

0.00 

0.00 

0.00 


0.00 

0.00 

0.00 

0.00 

Other (sim. only) 


0.16 

0.07 

3.12 


0.46 

0.09 

7 6 ? 

5.2! 

Total 

100 

100 

too 

1 00 

100 

100 

100 

100 

100 


B. Simulant Particle Size Distribution 

We used multiple sources of simulant size distribution data, and in most cases multiple data sources are 
represented per simulant. 

1 . ihy Sieving 

Some data are from dry sieving methods and reported by weight%. The data for OB-1 comes from Trow 
Analytical, Ltd. The analyses for JSC- 1 A and NU-LHT- 1M were performed in the lab of Dr. Susan Batiste at the 
University of Colorado Boulder. 

2. SEM and image processing 

We have size data from QEMSCAN*' SEM/EDS analysis, reported by weight%, for all simulants except for NU- 
LHT- 1 D. It should be said that grain mounts used for SEM imaging are polished and thus provide a sectioned 
sample, and that most particles will not be sectioned at their plane of greatest diameter equivalent. For this reason 
such results are sometimes referred to as an SSD (sectional size distribution) rather than a PSD. The high number of 
particles counted partially offsets this effect, but there will always be a slight bias towards finer particles in an SSD. 
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This can be partially compensated for by stereological techniques and we are pursuing this approach. For now. we 
caution the users to take this into account, but also remind them that all simulants were measured by this method and 
thus any problems will be consistent across that portion of the dataset. 

3. Liquid dispersion and laser dijfractometry 

We have data for NU-LHT-IM, -2M, and -ID, and JSC- 1 A from liquid dispersed laser diffractometrv. Susan 
Batiste at the University of Colorado Boulder measured NU-LHT-IM and JSC-1A, while the Bureau of Mines 
analyzed NU-LHT-2M and - ID. 

These data are presented as volume% rather than as weight%. If the particle composition distribution were 
consistent across the size fractions then the data would be equivalent, but this is not true for lunar regolith and is 
likely not to be true for simulants. However, we judge it likely that the deviations in density across the size fractions 
are of small effect. We leave it to the user to evaluate these ratings until more data are gathered and analyses are 
presented. Again, the method is consistent for the four simulants measured and thus is of comparative value. 

This analytical method yields more bins of data (smaller size fractions) than the FoM software allows. We have 
summed the bins to best match the bins in the literature lunar data 8 . 


V. Results of Figure of Merit Analyses 
A. Particle Type Composition Results 

All composition FoMs were run using a revised form of Figure of Merit Version I software not yet released. 
Table 2 presents the Figure of Merit composition results for all simulants tested against our 64001/64002 lunar 
reference material. 


Table 2. Results of Figure of Merit 
composition analysis. Figure of Merit V. 
1, Rev. 1 algorithm used with lunar 
reference material 6400 1 /64002. 

64001/64002 

simulant reference 


NU-LHT-1 M 

0.65 

NU-LHT-2M 

0.55 

OB-1 

0.28 

JSC-1 

0.33 

JSC-1 A 

0.35 

JSC-1 AF 

0.43 

MLS-1 

0.35 

FJS-1 

0.36 


8. Particle Size Distribution Results 

The particle size Figures of Merit were run using the Figure of Merit Version 1 software. The results are shown 
in Table 3. The range of size bins for 6400 1 '64002 ' is broader than for any of the simulants. All simulant PSD’s are 
compared to the entire 64001/64002 PSD and the results are shown in the first column of Table 3. Most of the 
simulant PSD’s only extend to -l mm. The second column of Table 3 shows comparisons of all simulants to the 
<!mm fractions of 64001/64002. These fractions were recalculated to sum to 100 weight%. For the two simulants 
specifically intended to be dust simulants, another normalized subset of literature data 8 was used for comparison, 
this time recalculating the <90 pm fraction to sum to 100 weieht%. We show the results for this subset in column 3 
of Table 3. 

it is a complex problem to consistently classify fragments of rock and breccia in lunar and terrestrial material, in 
lunar regolith, particles comprise a spectrum of varying glass content and fracture. For this reason, ai! rock 
fragments and breccias in the simulants are classified as lithic fragments and compared to the abundance of all rock 
and breccia fragments in the regolith. 
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Table 3. Figure of Merit size results for all simulants against 64001/64002 lunar reference material. Simulant 
datasets were compared against the bulk average of 64001/64002, the <lmm subset of the data and the <90 
p ms u bset o ft he d ata; both reference subsets were recalcui ated to 100%. Analytical method i s in parentheses . 



64001/2 bulk 

i 

! 

l 

c z. 

to 

A 

64001/2 average to 


average 

nun average 

90 pm 

OB- 1 (section image analysis) 

0.23 

0.54 


NU-LHT-IM (section image analysis) 

0.23 

0.58 


NU-LHT-2M (section image analysis) 

0.17 

0.48 


JSC-I (section image analysis) 

0.22 

0.53 


JSC-i A (section image analysis) 

0.25 

0.56 


JSC-I AF (section image analysis) 

0.06 

0.23 

0.60 

M.LS-1 (section image analysis) 

0.20 

0.29 


FJS-I (section image analysis) 

0.26 

0.45 


OB-1 (dry sieve) 

0.59 



NU-LHT-IM (dry sieve) 

0.26 

0.75 


JSC-I A (dry sieve) 

0.35 

0.74 


NU-LHT-2M (laser diffractometry) 

0.29 

0.82 


NU-LHT- 1 D (laser diffractometry) 



0.54 

NU-LHT- 1 M (laser diffractometry) 

0.26 

0.64 


JSC-1 A (laser diffractometry) 

0.28 

0.74 



Agglutinates are a member of the particle spectrum including lithic and breccia fragments, but we interpret them 
to be sufficiently unique in their properties and abundance as to be worth differentiating. Furthermore, their 
characteristics as irregularly shaped, often vesicular particles composed of minerals in a glass matrix makes it 
possible for the iDiscover software user to distinguish them based on results from the automated beam technology. 

Because the lunar regolith reference 64001/64002 is composed of -32 moda!% agglutinates and 31 moda!% 
lithic fragments, simulants that do not approximate these abundances will score a low composition FoM score. They 
may still be appropriate simulants for many purposes, by virtue of their chemistry, shape, or size distribution. 
Conversely, a simulant with appropriate abundances of these particles may be inappropriate for some uses. 

VI. Conclusion 

These results are discussed in further detail in the Lunar Regolith Simulant User’s Guide released by the 
Marshall Space Flight Center simulant project group. Further data on lunar highlands composition and PSD is still 
being gathered and compiled. We encourage users to contact the authors at Marshall Space Flight Center for 
advisement as to simulant use. We predict that the Lunar Regolith Simulant User’s Guide will be updated at least 
annually, but new information is available constantly. These evaluations are ongoing, as is Figure of Merit 
development. Most importantly, simulant development is continuing. 
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